Cosmological bounds on masses of neutrinos and 
other thermal relics 



Steen Hannestad 

Physics Department, University of Southern Denmark, Campusvej 55, 
DK-5230 Odcnsc M, Denmark, E-mail: hannestad@fysik.sdu.dk 

Abstract. With the advent of precision data, cosmology has become an extremely 
powerful tool for probing particle physics. The prime example of this is the 
cosmological bound on light neutrino masses. Here I review the current status of 
cosmological neutrino mass bounds as well as the various uncertainties involved in 
deriving them. From WMAP, SDSS, and Lyman-a forest data an upper bound on the 
sum of neutrino masses of 0.65 eV (95% C.L.) can be derived with any assumptions 
about bias. I also present new limits on other light, thermally produced particles. 
For example, a hypothetical new Majorana fermion decoupling around the electroweak 
phase transition must have m < 5 eV. 
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1. Introduction 

Neutrinos are the among the most abundant particles in the universe. This means that 
they have a profound impact on many different aspects of cosmology, from the question 
of leptogenesis in the very early universe, over big bang nucleosynthesis, to late time 
structure formation. 

At late times (T^T EW ) neutrinos mainly influence cosmology because of their 
energy density and, even later, their mass. 

The absolute value of neutrino masses are very difficult to measure experimentally. 
On the other hand, mass differences between neutrino mass eigenstates, (mi, 1112,1713), 
can be measured in neutrino oscillation experiments. 

The combination of all currently available data suggests two important mass 
differences in the neutrino mass hierarchy. The solar mass difference of <5m 2 2 — 8 x 1CT 5 
eV 2 and the atmospheric mass difference dm^ ~ 2.6 x 1CT 3 eV 2 [1-4]. 

In the simplest case where neutrino masses are hierarchical these results suggest 
that m 1 ~ 0, m 2 ~ 5m solar , and m 3 ~ 5m atmosphcric . If the hierarchy is inverted [5-10] 
one instead finds 1713 ~ 0, m 2 ~ <5m atm os P hcric, and m x ~ <5m atm ospheric- However, it is also 
possible that neutrino masses are degenerate [11-21], mi ~ m 2 ~ m 3 ^> 5m atm0 spheric, 
in which case oscillation experiments are not useful for determining the absolute mass 
scale. 

Experiments which rely on kinematical effects of the neutrino mass offer the 
strongest probe of this overall mass scale. Tritium decay measurements have been 
able to put an upper limit on the electron neutrino mass of 2.3 eV (95% conf.) [22]. 
However, cosmology at present yields a much stronger limit which is also based on the 
kinematics of neutrino mass. 

Very interestingly there is also a claim of direct detection of neutrinoless double 
beta decay in the Heidelberg-Moscow experiment [23,24], corresponding to an effective 
neutrino mass in the 0.1 — 0.9 eV range. If this result is confirmed then it shows that 
neutrino masses are almost degenerate and well within reach of cosmological detection 
in the near future. 

Neutrinos are not the only possibility for stable eV-mass particles in the universe. 
There are numerous other candidates, such as axions and majorons, which might be 
present. As will be discussed later, the same cosmological mass bounds can be applied 
to any generic light particle which was once in thermal equilibrium. 

Here I focus mainly on the issue of cosmological mass bounds. Much more detailed 
reviews of neutrino cosmology can for instance be found in [25,26]. 

2. Neutrino Decoupling 

In the standard model neutrinos interact via weak interactions with e + and e~. In the 
absence of oscillations neutrino decoupling can be followed via the Boltzmann equation 
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for the single particle distribution function [27] 

Tt- Ev T v = Cc ^ (1) 

where C co ii represents all elastic and inelastic interactions. In the standard model all 
these interactions are 2 <-> 2 interactions in which case the collision integral for process 
i can be written 

1 r d 3 p 2 d 3 p 3 d 3 p 4 
colwUlj 2Si 7 2£ 2 (27r) 3 2£ 3 (27r) 3 2£4(27r) 3 

x (27r) 4 5 4 (p 1 +p 2 -p3+P4)A(A,/2,/3,/4)5|M|2 2 _ >34)i , (2) 

where S\M\\ 2 ^ U ^ is the spin-summed and averaged matrix element including the 
symmetry factor S = 1/2 if there are identical particles in initial or final states. The 
phase-space factor is A(/i, f 2 , / 3 , U) = fsh(l - /i)(l - f*) ~ fM^ ~ ~ h)- 

The matrix elements for all relevant processes can for instance be found in 
Ref. [28]. If Maxwell-Boltzmann statistics is used for all particles, and neutrinos are 
assumed to be in complete scattering equilibrium so that they can be represented by 
a single temperature, then the collision integral can be integrated to yield the average 
annihilation rate for a neutrino 

1 F \9l + 9\)T\ (3) 



7T 3 



where 



9l + 9r= sin4 °w + + sin 2 6 w f for v e 



sin 4 9 W + (-- + sin 2 6 W ) 2 for v^ T . (4) 
This rate can then be compared with the Hubble expansion rate 

tf = 1.66<?y 2 ^- (5) 

to find the decoupling temperature from the criterion H = T\ t=Td . From this one 
finds that T D (u e ) ~ 2.4 MeV, T D (v^ T ) ~ 3.7 MeV, when g* = 10.75, as is the case in 
the standard model. 

This means that neutrinos decouple at a temperature which is significantly higher 
than the electron mass. When e + e~ annihilation occurs around T ~ rn e /3, the 
neutrino temperature is unaffected whereas the photon temperature is heated by a 
factor (11/4) 1 / 3 . The relation T^/T 7 = (4/11) 1 / 3 ~ 0.71 holds to a precision of roughly 
one percent. The main correction comes from a slight heating of neutrinos by e + e~ 
annihilation, as well as finite temperature QED effects on the photon propagator [28- 
41]. 

3. Neutrinos in structure formation 

Given that the neutrino temperature is known, the present contribution to the matter 
density from massive neutrinos is 

^ - "to- < 6 > 
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Thus, even a sub-eV neutrino mass gives a significant neutrino contribution to the energy 
density and therefore has an effect on structure formation. 

Perturbations in the neutrino distribution can be followed by solving the Boltzmann 
equation [42] through the epoch of structure formation. There are publicly available 
codes, such as CMBFAST [43], which can do this. 

The main difference between neutrinos and cold dark matter is that neutrinos are 
light and only become non-relativistic around the epoch of recombination, significantly 
later than matter-radiation equality. While they are relativistic, neutrinos free-stream 
a distance roughly given by A = 2n/k = cr, where r is conformal time. The 
total free-streaming length is therefore roughly given by A/ ~ ct(T = m v ). At 
scales smaller than this, the solution to the Boltzmann equation is exponentially 
damped. For scales larger than the free-streaming length neutrino perturbations are 
unaffected. It is therefore possible to divide the solution into two distinct regimes: 1) 
k > r(T = to): Neutrino perturbations are exponentially damped 2) k < r(T = to): 
Neutrino perturbations follow the CDM perturbations. Calculating the free streaming 
wavenumber in a flat CDM cosmology leads to the simple numerical relation (applicable 
only for T cq > to > T ) 



20 Mpc (T x 

Afs 



1 + log 3.9 




(7) 



In Fig. 1 transfer functions for various different neutrino masses in a flat ACDM 
universe (Q m + Q u + Q\ = 1) are plotted. The parameters used were fi& = 0.04, 
n CD M = 0.26 - Q v , tt A = 0.7, h = 0.7, and n = 1. 

When measuring fluctuations it is customary to use the power spectrum, P(k,r), 
defined as 

P(k,r) = \5 k (r)\ 2 . (8) 

The power spectrum can be decomposed into a primordial part, Po{k), and a transfer 
function T(k, r), 

P(k,T) = P (k)T(k,T). (9) 

The transfer function at a particular time is found by solving the Boltzmann equation 
for S(t). 

At scales much smaller than the free-streaming scale the present matter power 
spectrum is suppressed roughly by the factor [44] 

AP(k) AT(k,r = T ) o O, 



P(k) T(k,r = r ) ~n m ' (10) 

as long as VL V Q m . The numerical factor 8 is derived from a numerical solution of the 
Boltzmann equation, but the general structure of the equation is simple to understand. 
At scales smaller than the free-streaming scale the neutrino perturbations are washed 
out completely, leaving only perturbations in the non-relativistic matter (CDM and 
baryons). Therefore the relative suppression of power is proportional to the ratio of 
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Figure 1. The transfer function T(k,t = to) for various different neutrino masses. 
The solid (black) line is for m„ = 0, the long-dashed for m v = 0.3 eV, and the dashed 
for vn v = 1 eV. 

neutrino energy density to the overall matter density. Clearly the above relation only 
applies when VL V <C fl m , when fl u becomes dominant the spectrum suppression becomes 
exponential as in the pure hot dark matter model. This effect is shown for different 
neutrino masses in Fig. 1. 

The effect of massive neutrinos on structure formation only applies to the scales 
below the free-streaming length. For neutrinos with masses of several eV the free- 
streaming scale is smaller than the scales which can be probed using present CMB data 
and therefore the power spectrum suppression can be seen only in large scale structure 
data. On the other hand, neutrinos of sub-eV mass behave almost like a relativistic 
neutrino species for CMB considerations. The main effect of a small neutrino mass on 
the CMB is that it leads to an enhanced early ISW effect. The reason is that the ratio of 
radiation to matter at recombination becomes larger because a sub-eV neutrino is still 
relativistic or semi-relativistic at recombination. With the WMAP data alone it is very 
difficult to constrain the neutrino mass, and to achieve a constraint which is competitive 
with current experimental bounds it is necessary to include LSS data from 2dF or SDSS. 
When this is done the bound becomes very strong, somewhere in the range of 1 eV for 
the sum of neutrino masses, depending on assumptions about priors. This bound can be 
strengthened even further by including data from the Lyman-a forest and assumptions 
about bias. In this case the bound on the sum of neutrino masses becomes as low as 
0.4-0.6 eV. 
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Figure 2. A% 2 as a function of the the sum of neutrino masses m u . 

In Fig. 2 Ax 2 is shown for an analysis which includes the WMAP CMB data [45], 
the SDSS galaxy survey data [46,47], the Riess et al. SNI-a "gold" sample [49], and 
the Lyman-a forest data from Croft et al. [50] . For the Lyman-a forest data, the error 
bars on the last three data points have been increased in the same fashion as was done 
by the WMAP collaboration [45], in order to make them compatible with the analysis 
of Gnedin and Hamilton [51]. In addition to the neutrino mass, which I take to be 
distributed in three degenerate species, I take the minimum standard model with 6 
parameters: Q m , the matter density, Qf,, the baryon density, H , the Hubble parameter, 
and r, the optical depth to reionization. The normalization of both CMB, LSS, and 
Ly-a spectra are taken to be free and unrelated parameters. The priors used are given 
in Table 1. 

Table 1. Priors on cosmological parameters used in the likelihood analysis. 



Parameter 


Prior 


Distribution 


f2 = Q m -\- Vlx 


1 


Fixed 


h 


0.72 ±0.08 


Gaussian [48] 


n b h 2 


0.014-0.040 


Top hat 


n s 


0.6-1.4 


Top hat 


T 


0-1 


Top hat 


Q 




Free 


b 




Free 
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Table 2. Various recent limits on the neutrino mass from cosmology and the data 
sets used in deriving them. 1: WMAP data, 2: Other CMB data, 3: 2dF data, 4: 
Constraint on as (different in 4 a , 4 b , and 4 C ), 5: SDSS data, 6: Constraint on H , 7: 
Constraint from Lyman-a forest. 



Ref. 


Bound on J2 m u 


Data used 


Spergel et al. (WMAP) [45] 


0.69 eV 


l,2,3,4 a ,6, 7 


Hannestad [52] 


1.01 eV 


1,2,3,6 


Allen, Smith and Bridle [53] 


0.56lg;i 6 eV 


l,2,3,4 b ,6 


Tegmark et al. (SDSS) [47] 


1.8 eV 


1,5 


Barger et al. [54] 


0.75 eV 


1,2,3,5,6 


Crotty, Lesgourgues and Pastor [55] 


1.0 (0.6) eV 


1,2,3,5 (6) 


Seljak et al. [56] 


0.42 eV 


1,2,4 C ,5,6,7 


The present work 


0.65 eV 


1,5,6,7 



In this particular analysis, the 95% C.L. upper bound on the sum of neutrino masses 
is 0.55 eV. 

In Table 2 the present upper bound on the neutrino mass from various analyses is 
quoted, as well as the assumptions going into the derivation. As can be gauged from 
this table, a fairly robust bound on the sum of neutrino masses is at present somewhere 
around 0.5-1 eV, depending on the specific priors and data sets used. 

4. General thermal relics 

While light neutrinos are the canonical light, thermally produced particles, there are 
other species which can have the same properties. At T ~ M P1 gravitons are presumably 
in thermal equilibrium. However, since inflation occurs at a lower temperature, there 
should be no thermal graviton background present. 

Other examples include majorons, axions, etc. Any such species which was once in 
equilibrium is characterized by only two quantities, its mass, mx, and the temperature 
at which it decoupled from thermal equilibrium, Tr> At this temperature the number of 
degrees of freedom was g*. The relevant masses are <C MeV so that the particles decouple 
when they are relativistic, i.e. at decoupling they are characterized by a Fermi-Dirac or 
Bose-Einstein distribution of temperature Td. 

In the present-day universe, the new particles will be non-relativistic and contribute 
a matter fraction 

Q h 2 = ]]}x9x_ 9*u x ( fermions ) 1 (bosons) (11) 
183 eV g* V ; 3 V ; K J 

where gx is the number of the particle's internal degrees of freedom while g* v is the 
effective number of thermal degrees of freedom when ordinary neutrinos freeze out with 
g w = 10.75 in the absence of new particles. 
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In the late epochs that are important for structure formation, the momentum 
distribution of the new particles is characterized by a thermal distribution with 
temperature Tx that is given by 



In Ref. [57] CMBFAST was modified to incorporate a such thermal relics, both 
fermionic and bosonic. 

Here I present an updated analysis for the fermionic case with gx = 2 (i.e. a single 
Majorana fermion), using the same data sets as described above. Fig. 3 shows the 68% 
and 95% likelihood contours for VL x h 2 and g*. Overlayed are isocontours for particle 
masses. 

As an example, a species freezing out at the electroweak transition temperature 
has g* = 106.75, and from the figure it can be seen that the upper bound on the mass 
of such a particle is around 5 eV. For a corresponding scalar the mass bound will be 
slightly higher. For a species decoupling after the QCD phase transition where <^ 20 
the mass bound is roughly m <^ 1 eV. 




(12) 
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Figure 3. The full lines show 68% and 95% confidence regions in the (£l v h 2 , 
for the case where the thermal relic is a fermion with g = 2. 



<7») plane 



5. Conclusion 



I have reviewed the present status of cosmological mass bounds on neutrinos and other 
light, thermally produced particles. Already now these bounds are about an order of 
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magnitude stronger than current laboratory limits on the neutrino mass, albeit more 
model dependent. In principle the neutrino mass bound can be evaded either partially 
or completely by introducing broken scale-invariance in the primordial power spectrum 
[58] or by coupling neutrino strongly to a massless scalar so that they decay while still 
semi-relativistic [59]. 

In the coming years a wealth of new cosmological data will become available from 
such experiments as the Planck Surveyor CMB satellite. This is likely to allow for a 
measurement of the (sum of) neutrino masses in the 0.1 eV regime [60-63] . Together 
with direct detection experiments like KATRIN [64] and future neutrinoless double beta 
decay experiments, cosmology will provide the answer to whether neutrino masses are 
hierarchical. 
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